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The streamwise lift distribution of a wing-canard-stabilator-body configuration was 
varied to study its effect on the near-field sonic boom signature. The investigation was 
carried out via solving the three-dimensional Euler equation with the OVERFLOW-2 flow 
solver. The computational meshes were created using the Chimera overset grid topology. 
The lift distribution was varied by first deflecting the canard then trimming the aircraft with 
the wing and the stabilator while maintaining constant lift coefficient of 0.05. A validation 
study using experimental results was also performed to determine required grid resolution 
and appropriate numerical scheme. A wide range of streamwise lift distribution was 
simulated. The result shows that the longitudinal wave propagation speed can be controlled 
through lift distribution thus controlling the shock coalescence. 


Nomenclature 


c 

= 

wing root chord length 

Cd, C l C m 

= 

non-dimensional drag, lift, and pitching moment coefficients 

Cf 

= 

flat plate skin friction coefficient 

CGT 

= 

Chimera Grid Tools 

h 

= 

height (positive = above aircraft, negative = below aircraft) 

HLLC 

= 

Harten-Lax-Van Leer with Contact discontinuities 

/ 

= 

length of aircraft 

M 

= 

Mach number 

P 

= 

static pressure 

RANS 

= 

Reynolds Averaged Navier-Stokes 

Re 

= 

Reynolds number 

U 

= 

freestream velocity 

V 

= 

longitudinal cartesian coordinate 

y + 

= 

wall- spacing 

a 

= 

angle-of-attack (deg) 


= 

canard deflection angle (+ trailing-edge down) 

6 t 

= 

stabilator deflection angle (+ trailing-edge down) 

8 W 

= 

wing deflection angle (+ trailing-edge down) 

M 

= 

Freestream Mach angle (deg) 

Subscript 

ref = 

freestream reference value 


I. Introduction 

I NVESTIGATION of sonic boom has been one of the major areas of study in aeronautics due to the benefits a 
low-boom aircraft has in both civilian and military applications. Current Federal Aviation Administration 
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regulations prohibit supersonic flight over land due to potential effects the sonic boom may have on structures and 
humans. Consequently, numerous studies have been performed to investigate the characteristics of shock systems. 
The topic of special interest is the effect of streamwise lift distribution has on sonic boom. It is known that the lift 
affects the strength of the shock generated on the lifting surface, thus affecting the relative longitudinal propagation 
speed of the shock system. As the shock coalescence effect is the result of shocks having differential longitudinal 
propagation speed with respect to with other, it is possible to control this phenomenon via streamwise lift 
distribution. 

The effect of the streamwise lift distribution on the sonic boom was first investigated by Ferri 1 in 1969. Using a 
series of simple tandem wing configurations with varying distance and size, Ferri showed that the strength of the 
bow shock can be reduced by increasing the amount of lift carried near the nose of the configuration. This research 
was extended by Fomin 2 who investigated the influence of the streamwise lift distribution on the acoustic 
parameters using an aircraft model in a tandem wing configuration with different wing areas. The result showed that 
acoustic shock level can be reduced through proper distribution of lift along the length of the airframe. 

The effect of varying lift distribution on shock position and strength was investigated in flight in the Lift and 
Nozzle Change Effect on Tail Shock 3 (LaNCETS) project initiated at NASA Dryden Flight Research Center. Two 
aircraft were used for the project: 1) NASA aircraft NF-15B, tail No. 837 (NASA837) and 2) NASA aircraft F-15B, 
tail No. 836 (NASA836). NASA837, equipped with canard and variable area ratio nozzle, was utilized as the 
research aircraft while NASA836 served as the probing aircraft in the near-field. The lift distribution of the research 
aircraft was varied by deflecting the canard while maintaining trimmed level flight. 

The study presented in this paper directly supports the LaNCETS project with the computational analysis on the 
effect of streamwise lift distribution variation on shock. A simplified wing-canard-stabilator-body model was 
utilized instead of a full NASA837 in order to validate the computational tools, simplify the problem, and isolate the 
effects that may influence the sonic boom. The lift distribution was varied via canard deflection as it was done 
during flight. The angle of attack and the lift coefficient were kept constant for all configurations. In addition, all 
configurations were trimmed to steady, leveled flight. 

The computational approach is first described in this paper. A code validation study using a delta wing-body 
wind tunnel test case is described in Section III. The results of the lift distribution on the sonic boom signature are 
discussed in Section IV. 


II. Methods 

A. Flow Solver 

OVERFLOW-2 4 was chosen as the flow solver. OVERFLOW-2 is a viscous, compressible, three-dimensional 
Reynolds Averaged Navier-Stokes (RANS) solver specifically designed to utilize the structured Chimera overset 
grid 5 topology. It is a robust and comprehensive solver capable of employing various dissipation schemes, numerical 
schemes, flux limiters, and turbulence models. Although OVERFLOW-2 is a viscous flow solver, viscosity was not 
employed in majority of this study since viscous effects are negligible in sonic boom pressure signature 6 . However, 
viscosity was employed on conditions with numerical instability resulting from flow separations. The viscous 
calculations were done fully turbulent using wall-function of Spalart-Allmaras 7 turbulence model. The numerical 
scheme used throughout the study, the Harten-Lax-Van Leer with contact discontinuities 8,9 (HLLC) scheme, was 
chosen based on the scheme study performed in the validation section (Section III). Multigrid and grid sequencing 
algorithm were used to accelerate convergence rate. 

B. Computational Mesh 

The computational meshes were created using the Chimera overset grid 5 topology. This method simplifies the 
grid generation of geometrically complex multi -body configurations by allowing the user to construct the mesh from 
sets of relatively simple, overlapping, body-fitted grids. The process of grid generation was further streamlined by 
creating the mesh using the script library of Chimera Grid Tools 10 (CGT) v2.0. CGT is a comprehensive package 
containing variety of tools for creating grids using Chimera overset method. It contains tools for grid generation, 
manipulation, visualization, and diagnostic purpose. It can run in batch mode under the OVERGRID graphical 
interface or in scripting fashion using its extensive script library. The body-fitted computational meshes for viscous 
computation were created with non-dimensionalized wall-spacing, y + , of approximately +50. The y + was calculated 
using White’s 11 method shown in Equation 1 and Equation 2. 



Equation 1 


y' = R^^fy 

^ 0.455 

C r = — Equation 2 

f In 2 (0.06 Re x ) 

Re x denotes the Reynolds number at certain downstream distance, v (. x is typically set as the 10% of the reference 
length). The y denote the normal distance from the surface and C f denote the skin friction coefficient of a flat plate. 

The far-field grid was created aligned to the far-field Mach angle to improve the accuracy of the captured sonic 
boom. The far-field Mach angle was calculated using Eq. (3) 12 , shown below. The computational domain was 
extended out to 1 body-length in upstream and downstream direction and 2 body-lengths in spanwise direction. The 
domain connectivity was performed using Meakin’s X-Ray Object method 13 with double-fringe overlapping 
between all grids. 


/u = sin 



Equation 3 


C. Variation of Streamwise Lift Distribution 

A simple wing-canard-stabilator-body model was 
created to isolate the effect of lift distribution on 
shock. Propulsion devices such as inlet and nozzle 
were not modeled to limit the source of lift and shock 
systems. Features such as nose-boom, boat-tail, and 
wing-strakes were not modeled geometric simplicity. 

A symmetrical diamond airfoil with 5%c maximum 
thickness at 50% x/c was used for the wing, canard, 
and stabilator. All three lifting surfaces have 
approximately 45° leading-edge sweep. The fuselage Figure 1. Wing-canard-stabilator-body model 
was modeled as an axisymmetric body. The center of 

gravity was placed on the symmetry plane, 10.3 non-dimensionalized distances from the nose. The half-body model 
is shown in Figure 1. Some of the general features of NASA837 were conserved in the simplified wing-canard- 
stabilator-body model such as the span-to-chord ratio and the location of the individual lifting surfaces. 

As previously mentioned, the streamwise lift distribution was varied by deflecting the canard. Current method 
differs from that used by Ferri 1 and Fomin 2 which involved changing the size of the lifting surfaces to change the 
streamwise lift distribution. There are two advantages to changing the streamwise lift distribution via canard 
deflection. First, the aerodynamic center of individual lifting surfaces remains constant for different configurations. 
This eliminates effects that may result from having wings of different size. Second, the experiment is flight testable. 
Although it is possible to change the wing between the flight tests, the downtime and the cost associated the 
configuration change makes the test inefficient and impractical. The following procedure lists the steps for changing 
the streamwise lift distribution. 

1 . Deflect the canard 

2. Deflect the wing to reach the desired lift coefficient (C L = 0.05 ± 0.003) 

3. Deflect the stabilator to trim the configuration (C M = 0.0 + 0.005) 

4. Iterate between Step 2 to Step 3 until the configuration is trimmed at desired lift coefficient 

The deflections were considered positive for trailing-edge down. Lifting devices such as flaps and slats were not 
modeled for computational simplicity; instead, entire lifting surface was deflected as the control surface. A tool that 
proved to be very useful during this process was Missile DATCOM v.707 14,15 . Missile DATCOM is a preliminary 
design tool for estimating the aerodynamic forces and moment coefficient as well as control derivatives. It is capable 
of modeling missile-like geometries with various fins, nozzle, inlets, and protrusions. The lifting surface deflections 
were initially determined using Missile DATCOM v.707 then refined using OVERFLOW-2 inviscid results. 




The lift coefficient was kept constant at approximately 0.050 ± 0.003 for all configurations. The lift coefficient 
was chosen based on the requirement for numerical stability of the inviscid simulations. The configuration was 
considered trimmed if the moment coefficient was 0.000 + 0.005. The angle-of-attack was fixed at 1.5°. The altitude 
and Mach number was matched to one of the major test point flown in the LaNCETS project: Mach 1.4 at altitude 
40,000ft. 


III. Validation Study: 69° delta wing-body 

A numerical validation study was performed to 
determine the baseline grid resolution and investigate 
the performance of different numerical schemes. In 
addition, the force and moment coefficients from 
Missile DATCOM v.707 were compared to the 
inviscid result of OVERFLOW-2 for various angle-of- 
attacks. 

The “model 4” geometry described in a 1973 
wind-tunnel study by Huntun, Hicks, and Mendoza 16 
was selected as the validation model due to ample 
wind-tunnel data provided in the study. The model is 
comprised of an axisymmetric fuselage with a 69° 
leading-edge swept delta wing, shown in Figure 2. The 
airfoil of the wing is the symmetrical diamond airfoil with 5%c maximum thickness at 50% x/c. 

Due to lack of a sting geometry description, it was approximated by digitizing and extrapolating the data from a 
sketch of a wind-tunnel apparatus provided in the wind-tunnel study. 

The validation was performed against a non-lifting case (C L = 0.0, a = 0.0°) and a lifting case (C L = 0.08, a = 
2.56°) at freestream Mach number of 1.68. The freestream Mach angle, calculated using Eq.l, was 36.5° for the 
non-lifting case. The freestream Mach angle was adjusted for the lifting -case. The near-field pressure signature was 
measured parallel to the freestream velocity at h/l = +3.6 on the symmetry plane. The computational surface mesh 
and the volume mesh are shown in Figure 3. 




a) b) 

Figure 3. Validation model grid, (a) surface grid, (b) volume grid. 

A. Grid Resolution Study 

Grids of different densities were simulated to determine the base resolution required to accurately capture the 
propagating shock. The grid dimensions are tabulated in Table 1, shown below. 



Figure 2. “Model 4” wing-body validation model 16 . 


Grid 

Mach-Aligned 
Far-Field (JxKxL) 

Body (JxKxL) 

Wing (JxKxL) 

Sting (JxKxL) 

0.5x 

241x204x89 

110x25x25 

97x20x25 

60x61x25 

lx 

369x304x89 

193x50x50 

197x40x50 

136x61x50 

2x 

625x504x145 

373x 100x100 

397x80x100 

263x61x100 


Table 1. Grid Resolution for Validation Study 



The resolution of the Mach-aligned far-field grid was 
varied only in the regions of interest: the upstream and 
the region of shock propagation. Although this method 
would have been invalid for subsonic and transonic 
flows, it is valid in supersonic flow regime since the 
flow generally propagates toward downstream. It has 
been shown in a wind tunnel test that supersonic flow 
can propagate in the upstream direction within the 
boundary layer region 17 . However, for this test, it can 
be assumed that flow does not propagate upstream. 

The normalized propagated pressure of the 0.5x, 
lx, and 2x grids were compared against the wind- 
tunnel data in Figure 4. As shown in the figure, the 
effect from increasing grid resolution is negligible 
except in the body-sting connection region. However, 
the difference in the result between different grid 
resolutions is negligible. Overall the computational 
results are in agreement with the experimental data 16 
for a zero-lift condition. 

Comparison between the inviscid OVERFLOW-2 
simulation and the wind tunnel measurement for C L of 
0.08 is shown in Figure 5. The bow shock and the wing 
shocks show excellent agreement. An extra shock, 
shown in CFD but not present in wind-tunnel data, at 
1.1 < x/1 < 1.3 is the compression and the expansion of 
pressure between the wing and the tail. The cause of its 
existence in CFD may be due to inaccurate sting 
geometry or lack of viscous effects. Similar result for 
the tail is shown in other validation attempts done by 
Casper 18 and Wintzer 19 . Both Casper and Wintzer used 
CART3D compressible inviscid flow solver with grid 
adaptation. 

B. Numerical Scheme Study 

Different numerical schemes were compared to 
investigate the affects the schemes may have in 
capturing the shock. The non-lifting case was used for 
this portion of the study. The lx grid, described in 
Table 1, was chosen as the grid resolution of choice. 
The numerical schemes compared were Roes upwind, 
HLLC upwind 8 ’ 9 , AUSM+ 20 , and Yee’s Total Variation 
Diminishing (TVD) scheme 21 , shown in Figure 6. The 
numerical simulations were done using 3 rd order spatial 
accuracy for Roes upwind, HLLC upwind, and 
ASUM+. The TVD scheme was simulated using 2 nd 
order spatial accuracy. The HLLC scheme was also 
compared against the 5 th order Weighted Essentially 
Non-Oscillatory 22 (WEN05) and Weighted Essentially 
Non-Oscillatory Modified 22 (WEN05M) schemes, 
shown in Figure 7. 

As shown in Figure 6, Yee’s TVD has higher 
dissipation compared to Roes and HLLC scheme. It is 
also noticeable that ROES and HLLC scheme produce 
almost identical results. The 3 rd order AUSM+ scheme 
and the 5 th order WEN05 and WEN05M schemes 
produces unphysical undershoots and overshoots 



Figure 4. Validation grid resolution study result 
(C L =0). Experimental data from Ref. 16. 



Figure 5. Validation grid resolution study result 
(C l = 0.08). Experimental data from Ref. 16. 
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Figure 6. 2 and 3 order numerical 
comparison. Experimental data from Ref. 16: 
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Figure 7. 3 order and 5 order numerical scheme 
comparison. Experimental data from Ref. 16: 


around the shock. These effects may be partially mitigated by doing the full viscous simulation. 

The scheme chosen for the rest of the study, based on the scheme comparison results, is the HLLC scheme. 
Although HLLC and Roes upwind scheme produced similar results, HLLC was chosen due to higher theoretical 
robustness compared to the ROES upwind scheme 8 . 

C. Missile DATCOM v.707 and OVERFLOW-2 Comparison 

The lift, drag, and pitching moment coefficients of 
the validation model estimated using Missile 
DATCOM v.707 were compared against that from 
OVERFLOW-2 inviscid simulation. The lx grid, 
described in Table 1, was utilized for the 
OVERFLOW-2 simulation. The comparison was 
necessary since Missile DATCOM was used to 
initially determine the lifting surface deflection of the 
model used in the streamwise lift distribution study. 

The lift and pitching moment coefficients were 
compared for various angles-of-attack. The result is 
plotted in Figure 8. The result shows excellent 
agreement between Missile DATCOM v.707 and 
OVERFLOW-2 for low angle-of-attack. Although the 
range of angle-of-attacks compared is limited, it is 
sufficient for this study. 



Figure 8. Missile DATCOM and Overflow-2 inviscid 
result comparison: C L , C D , C M . 


IV. Streamwise Lift Distribution Variation Study: wing-canard-stabilator-body 

Total of four configurations were analyzed in this study. The configurations are summarized in Table 2, shown 
below. Configurations 1, 2, and 3 were tested in flight; Configuration 4 was added to the test matrix for 
completeness and Configuration 0 was utilized only for the grid independence study. 


Configuration # 

Canard Deflection (deg) 

Wing Deflection (deg) 

Stabilator deflection 
(deg) 

0 

0.0 

0.0 

0.0 

1 

-5.5 

0.3 

-6.0 

2 

0 

-1.1 

0.0 

3 

+3.3 

-1.4 

1.0 

4 

+5.5 

-1.8 

3.0 


Table 2. Descriptions of test configurations. 


The lift generated by individual lifting surfaces for 
different configurations are shown in Figure 9. As 
previously mentioned, the lift coefficient was constant 
(C L = 0.05) and angle-of-attack was fixed at 1.5° for 
Configurations 1 through 4. Comparing the 
configurations, Configuration 1 generates largest 
positive lift with wing and largest negative lift with 
canard and stabilator. Configuration 2 generates 
positive lift on wing, canard, and stabilator. However, 
the lift on wing is significantly less than that generated 
by Configuration 1. Configurations 3 and 4 generate 
negative lift with wing and positive lift with canard 
and stabilator. It is clear to see that configurations 
tested span a wide, if not full, range of lift that could 
be distributed over an aircraft. 

The flow around Configuration 4 was resolved 
with viscosity due to flow separation at the stabilator. 
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Figure 9. Lift generated by individual lifting surfaces 
for different configurations. 


Configurations 1, 2, and 3 were computed without viscosity. All propagated pressures were measured parallel to the 
freestream velocity, at M=- 1.6 on the symmetry plane which is approximate location of the probing aircraft during 
the flight. The freestream Mach angle, at Mach 1.4, was 45.6°. 


A. Grid Independence Study 

The grid independence study was performed to 
determine the grid resolution required to accurately 
capture the propagated shock. Configuration 0 was used 
for this study. Due to a large number of grid points in 
the Mach-aligned far-field grid, the grid independence 
study was done separately for the body-fitted grids and 
the mach-aligned far-field grid. The body-fitted grid and 
the Mach-aligned far-field grid were tested at angle-of- 
attack of 1.5° and 2.0°, respectively. The grid 
independence was established based on the pressure 
instead of the force and moment coefficients since this 

study strictly focuses on the pressure signature. The grid 
dimensions for the body-fitted grids and the Mach- 
aligned far-field grids are tabulated in Table 3 and Table 
4, respectively. The grid independence of the body- 
fitted grids was established by comparing the surface 
pressure of 0.5x, lx, and 2x grids, as shown in Figure 
10. It is clear from the figure that surface pressure does 
not change with the increase in the grid resolution, thus 
establishing grid independence at 0.5x grid resolution. 
Similarly, the grid independence of the mach-aligned 
far-field grid was established by comparing the 
propagated pressure for lx, 2x, and 3x grid, shown in 
Figure 11. The grid independence was established at lx 
grid resolution both body-fitted and far-field. 



Figure 10. Wing-canard-stabilator-body body-fitted 
grid independence study. 



Figure 11. Wing-canard-stabilator-body body-fitted 
grid independence study. 


Grid 

Bodys 

(JxKxL) 

Lifting Surfaces 
(JxKxL) 

Collar 

(JxKxL) 

0.5x 

241x30x25 

97x20x25 

97x19x35 

lx 

280x50x50 

197x40x50 

197x39x50 

2x 

348x 120x50 

397x159x100 

397x79x100 


Table 3: Body-fitted grid study dimensions 


Grid 

Far-field (JxKxL) 

lx 

297x225x105 

2x 

553x455x 105 

3x 

809x655x105 


Table 4: Mach-aligned far-field grid study dimensions 

B. Effect of Streamwise Lift Distribution Variation on Shock 

The normalized propagated shock is plotted in Figure 12 for different streamwise lift distributions represented 
by various canard deflections. The lift generated by individual lifting surfaces for different configurations are shown 
in Figure 9. More details of shock system are visible for Configuration 4 due to the effect of viscosity. 

As mentioned briefly in the Introduction section, the compression and expansion waves of a large shock system 
propagate faster in forward and aft direction, respectively, than that of a smaller shock system. Consequently, the 
phenomenon of shock coalescence occurs because the waves of a larger shock system catches up to that of a smaller 
shock system. In addition, it is known that the strength of the shock generated increases with the increase in lift 
generated by the lifting surface. Thus it is possible to control relative speed of the shock propagation through 
engineered distribution of streamwise lift. Consequently, it is possible to control the shock coalescence. 


A conventional wing-tail configuration aircraft 
generates most of the lift with the wing while using the 
tail as the trimmer. Thus the shock generated on the 
wing is significantly stronger than that generated on the 
tail. As shock systems propagate toward the ground, 
the compression wave from the wing catches up to the 
nose shock system and coalesce with it. Similarly, the 
expansion wave from the wing coalesces with the 
shock system from the tail. Consequently, although 
many shock systems were generated on the aircraft, 
only two shocks are heard on the ground; the bow 
shock and tail shock. A similar shock system is shown 
in Figure 12 on Configuration 1. 

Examining the lift distribution and the pressure 
signature of Configuration 1, the wing generates large 



Figure 12, 
stabilator-body 
configurations. 


Propagated sonic boom of wing-canard- 
at h/l = -1.6 for different 


amount of positive lift while the canard and stabilator generates negative lift. Such lift distribution generates a large 
shock system which propagates much faster than the shocks system generated on the canard and the tail. The figure 
shows that difference in speed of propagation is large enough for the compression wave of the wing to coalesce with 
the shock system of the canard within 1.6 body-lengths. It is also noticeable that expansion wave from the tail has 
coalesced with the expansion wave of the tail. 

However, it is also possible to prevent shocks from coalescing. Knowing that shock propagation in forward and 
aft direction is proportional to the strength of the shock, one can prevent shock coalescence by generating shocks of 
approximately same strength for all lifting surfaces. Such shock system would not coalesce while propagating 
toward the ground since all shock systems would be propagating forward and aft at same speed. Such case is shown 
in Figure 12 on Configuration 3. 

Configuration 3 has its wing generating negative lift while the canard and the stabilator are generating positive 
lift. This lift distribution generates shock systems of approximately equal strength for the wing, canard, and the 
stabilator. It is shown in figure that shocks of the three lifting surfaces have not coalesced within 1.6 body-lengths. 
In addition, the expansion wave from the tail is still a separate shock from the expansion wave of the tail. It is also 
noticeable that strengths of shock system of wing, stabilator, and canard are still approximately equal. As the shock 
systems propagate toward the ground, the shock system of the canard is likely to coalesce with the shock system of 
the nose. Similarly, the shock system of the stabilator is likely to coalesce with the shock system of the tail. 
However, unlike the conventional aircraft, three relatively weaker shock systems will hit the ground instead of a 
single strong shock system. Configurations 2 and 4 show intermediate shock system compared to Configurations 1 
to 3. Same trends were shown in the LaNCETS flight data. 

It can be concluded from above observation that it is possible to prevent shocks from coalescing with one 
another. By default the wing generates the most amount of lift thus produce strong shock compared to canard and 
stabilator. Thus as the shock systems propagate from the aircraft toward the ground, the shock system of individual 
lifting surfaces coalesce due to different forward and aft propagation speed. However, with careful distribution of 
lift, it is possible to various lifting surfaces to generate shocks of approximately equal strength which would have 
equal forward and aft propagation speed, thus preventing shock interaction. 


V. Conclusion 

The effect of change in streamwise lift distribution on sonic boom and the shock structure was investigated 
using a simple wing-canard-stabilator-body configuration. The study demonstrated the ability to control the shock 
coalescence by controlling the shock system’s relative longitudinal propagation speed using streamwise lift 
distribution. The result also showed that it could be more advantageous for aircraft to offload the lift from the wing 
to stabilator and canard to prevent shocks from coalescing in the near-field region. An extensive validation study 
was performed to determine the base grid resolution and the effects of different numerical schemes which increase 
the confidence in the result. 
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• Lift and Nozzle Change Effect on Tail Shock 

- Engine plume effects (AIAA-2009-1054) 

- Streamwise lift distribution effects on shock coalescence 



Method: Streamwise Lift Distribution Variation 


• Simple wing-canard-stabilator-body configuration 

• Missile DATCOM v.707 for lifting surface 
deflections 

• OVERFLOW-2 for checking and refining 

• Inviscid simulation unless viscosity needed 


NO 


Deflect canard 


Deflect wing and 
achieve desired lift 
coefficient 



YES 


Done 


Tools 
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OVERFLOW-2 Flow Solver 

• 3-Dimensional, compressible, 
viscous, Reynolds Averaged 
Navier-Stokes Solver 

• Overset/chimera grid topology 

• Multiple turbulence models 

— wall-function variation 

- DES implementation on some 
turbulence models 

• Steady-State and time- 
accurate simulation 

• Convergence acceleration 
schemes 

• 6-DOF simulation 


Missile DATCOM v.707 

• aerodynamic prediction 
code 

• Capable of modeling 
missile-like geometries with 
multiple fin stations, 
protrusions, nozzles, and 
noses 

• Wide range of flight regime 

• Validated against wind 
tunnel as well as CFD 


Validation: 
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Description 
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69° delta wing-body 

- Symmetrical 5%c max. 
thickness at mid. chord 
diamond airfoil 

- wind tunnel data provided 
NASA TN D-7160 



- C L = 0.00 (0°) validation 

- C L = 0.08 (2.56°) validation 

- Mach 1.68 

Tests: 

- Grid Resolution Study 

- Numerical Scheme Study 

- Missile DATCOM / OVERFLOW-2 




comparison 


Validation: Grid Resolution Study 


• Zero-lift condition 


Grid 

Mach-aligned Far-Field 
(JxKxL) 

Body 

(JxKxL) 

Wing 

(JxKxL) 

Sting 

(JxKxL) 

Total # 
Grid Point 

0.5x 

241x204x89 

110x25x25 

97x20x25 

60x61x25 

~4.5M 

lx 

369x304x89 

193x50x50 

197x40x50 

136x61x50 

-11.3M 

2x 

625x504x145 

373x100x100 

397x80x100 

263x61x100 

-54.7M 


d/fd-d) 
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Validation: Grid Resolution Study Result (C L = 0) 


CD 
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Validation: Numerical Scheme Study 


Zero-lift condition 
Schemes tested: 

- 3 rd order Roes Upwind 

- 3 rd order Advection Upstream Splitting Method (AUSM+) 

- 3 rd order Harten-Lax-Van Leer with Contact discontinuities 
(HLLC) Upwind 

- 2 nd order Total Variation Diminishing (TVD) 

- 5 th order Weighted Essentially Non-Oscillatory (WEN05) 

- 5 th order Modified Weighted Essentially Non-Oscillatory 
(WEN05M) 
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000 
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015 
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dation: Numerical Scheme Study Result 
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Validation: Numerical Scheme Study Result 


<u 



CL 

CL 



d/Cd-d) 
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Validation: Lifting Case Result (C L = 0.08) 
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Validation: Missile DATCOM v.707 vs. OVERFLOW- 2 



2.0 


4.0 


6.0 


Validation: Result 

ROES upwind and HLLC upwind schemes produces nearly 
identical result 

Total Variation Diminishing scheme shows higher 
dissipation compared to other schemes 

AUSM+, WEN05, WEN05M produces too much 
unphysical oscillations to be useful for this test 

ROES upwind, HLLC upwind, AUSM+, WEN05, WEN05M 
agrees well with wind-tunnel measurement 

Missile DATCOM agree well with OVERFLOW-2 for 
small angle-of-attack 

HLLC chosen for this study 

lx grid resolution chosen 




Streamwise Lift Distribution: Description 

• Simple wing-canard-stabilator body model 


- Propulsion effects (inlets, nozzles) not modeled 

- Boat-tail, nose-boom, strake not modeled 

- Simplify and accelerate computation and grid generation 

- ~ 11 million grid points, 90% in the Mach-aligned far-field grid 

• Some of LaNCETS research A/C features conserved 

- Lifting surface span-to-chord ratio 

- Lifting surface position relative to nose 


17.62 


13.736 


CGx = 10.3 


8.62 



18.52 


Streamwise Lift Distribution: Test Matrix 


Configuration # 

Canard Deflection 
(deg) 

Wing Deflection 
(deg) 

Stabilator Deflection 
(deg) 

1 

-5.5 

0.3 

-6.0 

2 

0 

-1.1 

0.0 

3 

3.3 

-1.4 

1.0 

4 

5.5 

-1.8 

3.0 


*all deflections positive for trailing-edge down 


• Configuration 4 simulated with Spalart-Allmaras wall-function 
turbulence model due to flow separation (y + = 50) 

• Mach 1.4 

• Constant Angle-of-Attack = 1.5° 

• Constant lift coefficient = 0.05 

• Trimmed flight (zero pitching-moment coefficient) 


Lift Coefificient (CL) 


9 



ryden Flight Research Center 




Streamwise Lift Distribution: Result 


0.06 

0.05 

0.04 

0.03 

0.02 

0.01 

0.00 

- 0.01 

- 0.02 




■ Config 1 


■ Config 2 


■ Config 3 


■ Config 4 



-0.03 


Lifting Surfaces 


d/fd-d) 
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Streamwise Lift Distribution: Result 


0.050r 


0.025- 
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3 0.000 


-0.025 
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- 0.2 


canard 


— config 1: large pos. lift on wing 
— config 2: small pos. lift on wing 
— config 3: small neg. lift on wing 
— config 4: large neg. lift on wing 



0.2 0.4 0.6 0.8 1.0 1.2 

x/l [normalized distance from freestream Mach line] 


Conclusion 
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Detailed validations of OVERFLOW-2, Missile 
DATCOM v.707 presented 

The local Mach angle can be varied based on 

the streamwise lift distribution 

- Possible to prevent shock coalescence by 
distributing lift over the aircraft to create shocks of 
equal strength 


Questions? 
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